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We propose an implementation of external homogeneous magnetic fields in k-p Hamiltonians 
for holes in heterostructures, in which we made use of the minimal coupling prior to introduce 
the envelope function approximation. Illustrative calculations for holes in InGaAs quantum dot 
molecules show that the proposed Hamiltonian outperforms standard Luttinger model [Physical 
Review 102, 1030 (1956)] describing the experimentally observed magnetic response. The present 
implementation culminates our previous proposal [Phys. Rev. B 82, 155307 (2010)]. 

PACS numbers: 73.21.-b, 73. 21. La, 75.75.-c 
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I. INTRODUCTION 

The interaction of a magnetic field with a charged par- 
ticle with spin comes into the Hamiltonian through cou- 
pling to the total (orbital plus spin) angular momentum. 
In crystals, the total angular momentum is the sum of the 
Bloch angular momentum, which contains atomic orbital 
and spin contributions, and the envelope orbital angular 
momentum-*^ Determining the coupling constant (g fac- 
tor) is an important requirement to study the magnetic 
properties of materials. In bulk systems, the value of the 
g factor is strongly influenced by band mixing, spin-orbit 
interaction and crystal anisotropy^ 

In the last years, there is increasing interest in con- 
trolling and exploiting the g factor of carriers confined 
in semiconductor quantum dots (QDs) for spin prepa- 
ration, conservation and manipulation aiming at spin- 
tronic and quantum information devices.— ~— The mag- 
netism of these structures is significantly different from 
that of bulk crystals because the weak magnetic confine- 
ment is supplemented by a strong spatial confinement. 
The latter has a profound effect on the energy structure, 
carrier-carrier interactions, band mixing and spin-orbit 
interactions, 3 which, in turn, influence the g factor value. 
As a matter of fact, it has been shown that quantum con- 
finement in QDs leads to strongly anisotropic g factors for 
both electrons and holes fiir— as well as to a quenching 
of the g factor value J^r— 

Up to date, most theoretical studies investigating the 
magnetic response of QDs rely on effective mass and k • p 
models. The standard inclusion of magnetic fields in 
k • p Hamiltonians consists in replacing the canonical 
momentum p by the kinetic momentum p — qA, sup- 
plemented with the spin Zeeman term, in the differential 
equation fulfilled by the envelope function (here q is the 
carrier charge and A the potential vector defining the 
magnetic field) This approximation, hereafter re- 

ferred to as the Luttinger approximation, has been suc- 
cessful in explaining several experimental observations in 
heterostructures , 3 ' 20 ! 21 it is implemented in the widely 
employed semiconductor software package nextnano 22 ^ 23 
and it is currently being used to determine the g fac- 



tors of confined carriers ! 14 ' 23 ' 24 However, a number of 
situations have been identified where it provides quali- 
tatively incorrect predictions. For example, in quantum 
rings under axial magnetic fields, the Luttinger approx- 
imation predicts the optical gap to decrease with the 
field strength^ contrary to magneto-photoluminesence 
observations^ Similarly, in vertically coupled QDs, the 
Luttinger approximation predicts that an axial magnetic 
field can tune the hole tunneling rate^i but this effect 
is not observed in related experiments^ The underly- 
ing reason is that the Luttinger approximation includes 
off-diagonal magnetic terms in the Hamiltonian, which 
artificially enhance the heavy hole-light hole (HH-LH) 
band mixing i 25 ' 28 

In Ref. |29| . an alternative formulation of the mag- 
netic interaction was suggested, in which the replacement 
of the canonical momentum by the kinetic one is car- 
ried out prior to the envelope function approximation 
(EFA). The resulting Hamiltonian has no off-diagonal 
magnetic terms directly coupling HH and LH, and the 
results become then consistent with the experimental 
measurements i 25 ' 28 

In the present work, we extend Ref. theory in or- 
der to account for the spin Zeeman term, and iden- 
tify the coefficients that should accompany the mag- 
netic terms in this approximation, which were pending 
clarification. 28 The remote band influence is considered 
through the zero-field effective masses, which are known 
to provide a meaningful description even in strongly con- 
fined QDs2°-~— , and effective ^-factors. We run calcula- 
tions comparing this model with the Luttinger approxi- 
mation and show that the magnitude of the Zeeman split- 
ting we estimate is closer to experimental values for ver- 
tically coupled InGaAs QDs. Thus, our theory offers im- 
proved accuracy in current attempts to understand and 
predict g factor values of holes in QDs. 

The paper is organized as follows. In Section [ill we 
derive the multiband k-p-EFA Hamiltonian for holes in 
the presence of a magnetic field. Starting from a gen- 
eral formulation, the Hamiltonian is then particularized 
to the case of QDs with axial symmetry. In Section Hill 
we use the obtained Hamiltonian to calculate the Zee- 
man splitting of vertically coupled QD molecules. The 
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results are compared with previous implementations of 
the magnetic field and experimental data. Conclusions 
are given in Section IIVI 



where [7^, Kj] — iTiTTj — TTjiti. 

We may define er with components = , where k 
represent a cyclic permutations of x, y, x. Then, we have: 



II. THEORY 



The Hamiltonian 



The classical Hamiltonian of a charged particle with 
anisotropic mass response to external forces, subject to 
the action of a magnetic field defined by a potential vector 
A is: 



2 —(pi-qAi) 2 



2m,- 



(1) 



where q is the charge and m», i — x,y,z the anisotropic 
mass. 

Elemental particles, in addition to charge, have spin. We 
can introduce heuristically spin by making the formal 
replacement tt — > er • ir in the above equation and taking 
into account the next two identities involving vectorial 
operators: 



(cr ■ a) (it • b) = a ■ b + icr ■ (a x b) 



p x A + A x p 



-ihV x A 



-ihB 



(2) 



(3) 



In the above equations a, b are vector operators, the 
components of er are the Pauli matrices and B represents 
the magnetic field. This formal replacement turns the 



kinetic energy T 



_ (p-gAj 



into 



(p-gA) 2 qh 
1 D = er ■ ±3, 

2m 2m 



(4) 



as it should appear in the Dirac equation. If the mass is 
anisotropic we should write instead: 



'3"3 



(5) 



1 1 i 

— (er • 7r) (<t • 7r) = — 7f 2 + — er • 7F x 7F 



where 7r, 



The er • 7r x 7r term in (J5)) can be expanded as follows: 



1 



1 



CT • 7T X 7T = Or, 



Fj/,7T*J + 0-3/ 



^Jm y m 
1 . 



^m x m z 



(6) 



er-7rX7T = er-(7rx 7r). 



(7) 



In the other hand, since 7r x 7t = (p — qA) x (p — qA) — 
—q(p x A + A x p) = iqhV x A = iqh'B we find out 
that, 



% qh 

— er •7TX7T = er • B. 

2 2 



(8) 



From now on, we will restrict ourselves to the case of axial 
symmetry, i.e., to the particular case B = Bph, m x — 
m„ — mi. Then, the above term turns into — „ qh <j z Bq 

y - 1 - 1 2m± " u 

and the complete Hamiltonian reads, 



W = (Pa - qA a )-^— (p Q - qA a ) 



a— -L,z 

— — a z Bq = T — a z Bo 

2m± 2m±_ 



2m a 
qh 



(9) 



By employing the symmetric gauge A = ^ [— y, x,0], 
we see that A z = 0. Then, z component of the kinetic 
energy is not affected, T z — p z t^- P z > while in-plane 
component is: 

_ 1 q 2 A 2 , qA± q A ± . . 

7j_ = Pa_ P± H ■ P_l ■ (10) 

2m_i_ 2mj_ 2m_i_ 2 m^ 

Since m±(p,z) then, px(^j-) = ~[x, y, 0] ^(^) and 

p± ■ ■^ J ~ 5* = ■^ J =- • p±ty, so that the in-plane component 
of the Itinetic energy results 

1 q 2 A\ 



71 = Pa 



2m 



■PA 



2m 1 m 



A± ■ pa_ 



-n Vj_ V_l + — 7; {xpy - yp x ) 



2m 



8m± 2m± 



(11) 



and the complete Hamiltonian is in turn 

q 2 B 2 p 2 qBp > g^o 
8m_L 2mj_ 2 2m j_ 



a— _L,2 



(12) 

This equation particularized to effective electrons (m > 
0, q = -1) is, 

n e = - J2 V a J-V a +|^ + -^(4+0 (13) 
* — ' 2m Q ami 2mj_ 



CL — A^^Z 



What about holes (m = — |m| < 0, g = 1)? Holes are 
tricky particles that require a careful tackle. We know 
that in the case of a one-band model, electrons and holes 
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energy dispersion are mirror image of each other. Then, 
we should assume that, 



a=.L,z 



R2 „2 



2mJ 



Bp 
2\m, 



(L z + a z ) 
(14) 



B. Envelope and Bloch functions 

In solid state physics the wave function |W(r)) is 
expressed as a sum of products |W(r)) = \JJ z )i\f)i, 
where \JJ Z ) are the Bloch band-edge and |/) are the 
envelope functions Let us call H° to the Hamil- 
tonian, eq. (1141) . in absence of magnetic field and 
to the second and third terms of this equation 
describing the action of the external magnetic field. 
The Luttinger-Kohn H° matrix Hamiltonian operator 17 
acting on the envelope vector function can be obtained 
by applying H° onto |\&(r)). Then, left- multiplying by 
the different Bloch functions (JJ z \i and integrating over 
the unit cell. Afterwards, the effect of remote bands is 
incorporated by replacing the actual mass by effective 
masses in the matrix elements of H° . 

In the presence of magnetic field H° must be supple- 
mented by H^- 8 ' coming from and |\&(r)), through 
a similar procedure. Since all terms in THS B ^ act as pure 
multiplicative operators on the envelope function compo- 
nents except for L z , we have thatjSS 



n^\jj z )\f) = \m^\jj z 



\JJ Z 



Bp 
'2\m± 



L z \f) (15) 



the third term in becomes —k^lbBq^ z . The proce- 
dure yields the following non-zero matrix elements for the 
6x6 valence Hamiltonian for zinc-blende crystals, which 
includes heavy hole, light hole and split-off bands. 



4? 



4? 



L 44 



AB 
L (i(i 



-(71 +72) 
-(71 -72) 
-(71 -72) 



^ + f (F.-3/2) 
Blp 2 , Bo 



+ ^(F,-l/2) 



Blp 2 , Bo 



+ ^(F z + l/2) 



-(71 +72) 

'B 2 oP 2 , Bo 



Blp 2 , Bo 



+ ^(F + 3/2) 



-Kp B B 



- -Kp, B B 



+ -kUbBo 



= -71 



= -71 



8 +T(**.-l/2) 
^ + ^. + 1/2) 



1 ' R 

- -k p B B 

+ ^k'ubBq (18) 



where fig = \q\/2m is the Bohr magneton, while k and k' 
are effective g factors for holes which become 4/3 and 2/3 
respectively if we remove the contribution of the remote 
bands. 

Note that the magnetic terms in Eq. (fl8l) differ from 
those of the Luttinger approximation (see Table 3 in 
Ref. I29I ). In particular, there are no off-diagonal mag- 
netic terms mixing HH and LH subbands. They also 
differ from our previous proposal^ in two aspects: (i) 
the spin degree of freedom is now included, and (ii) the 
remote bands contribution to the linear-in-B term com- 
ing from the Bloch function (third term in Eq. (|17t ) is 
now included through hole g factors («, n). 



III. ILLUSTRATIVE CALCULATIONS 



Axially symmetric systems have well defined z- 
component F z of the total angular momentum and, addi- 
tionally, the components of the envelope function associ- 
ated to the Bloch function \JJ Z ) have also a well defined 
M = (F z — J z ) orbital angular momentum; 35 ' 36 so that 
L z \f) = (F z - J z )\f). Then, we calculate the (J' J' z , JJ Z ) 
matrix element of as follows: 



{j'j' z \n (B) (\JM)\f)) = 
(J'j' z \(n^ 



Fz J z 

2|m±l 



H< B > 



2 „2 



F z -J z 



8\mj_\ 2\mj_\ 



Bnl- 



B )\JJ z )-\f) (16) 



(L,+rT z ) (17) 



B n 



2 mi 



Next, we incorporate de effect of the remote bands like in 
electrons: the mass m± arising in the two first terms of 
H' B ' is replaced by the effective mass parameter appear- 
ing in the corresponding matrix elements of H°, while 



In this section we implement the magnetic terms de- 
scribed above in a 4-band k-p Hamiltonian coupling HH 
and LH states. Hereafter we refer to this as W K . For com- 
parison, we also implement the magnetic terms following 
the Luttinger approximation, Hi u tt (Table 3 of Ref. |23 
but adding diagonal spin terms, — k[XbBqJ z ), and our 
previous pro posal T-L m ass (four-band Hamiltonian taken 
from Ref. 12^. 

To test the performance of the different Hamiltonians, 
we compare with available experimental data for an In- 
GaAs/GaAs QD molecule subject to a longitudinal mag- 
netic fields The Hamiltonian is solved numerically for 
a structure formed by two vertically stacked cylindrical 
QDs. The QDs have radius R = 15 nm and height H = 2 
nm, with an interdot barrier of thickness S — 2.8 nm. 
InGaAs Luttinger parameters are used for the masses, 
(71 = 11.01, 72 = 4.18, and 73 = 4.84)^ and a valence 
band offset of 0.2 eV is taken at the interfaces. In bulk, 
the k constant takes a value of 7.68 (for pure InAs)^ 
However, quantum confinement severely quenches this 
value. In QDs, one can safely disregard the contribu- 
tion from remote bands, and simply consider that com- 
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ing from the HH-LH subband coupling^ Therefore, we 
take K = 4/3. 



> 
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FIG. 1: Magnetic field dispersion of the highest hole states 
in a QD molecule calculated with different implementations 
of the magnetic terms, (a): Hiutt, (b): H mass , (c): H K . 
Solid and dashed lines are used for the bonding and anti- 
bonding hole states. The arrow in (a) indicates the bonding- 
antibonding ground state reversal. 

Figure Q] shows the energy of the highest valence band 
states. These are the \Fz\ = 3/2 hole states with bond- 
ing (solid lines) and antibonding (dashed lines) molec- 
ular character^ Panels (a), (b) and (c) correspond to 
estimates obtained with Hiutt, Hmass, and H K , respec- 
tively. One can see there are conspicuous differences in 
the energy spectra. For example, both Hiutt and Hmass 
predict that for the ground state, the £?-linear term dom- 
inates over the _B-quadratic (diamagnetic) one, so that 
its energy increases with the field. This would imply a 
decrease of the excitonic gap, in sharp contrast with pho- 
toluminescence experiments of InGaAs QDs, where the 
gap increases quadratically, indicating that the diamag- 
netic term is dominant. This is precisely the situation 
predicted by H K , Fig.QJc). 




FIG. 2: Bonding-antibonding energy splitting as a function of 
the magnetic field. Dashed line: Hiutt- Dashed-dotted line: 
Hmass- Solid line: H K . 

Further insight is obtained by comparing energy dif- 
ferences within each spectrum. We first compare the 
energy splitting between the bonding and antibonding 



states, Abab as a function of the magnetic field. Fig. [2] 
shows A oa b calculated with the three Hamiltonians. Hiutt 
(dashed line) predicts that the energy splitting decreases 
with Bq, becomes zero at Bq = 5.3 T and negative after- 
wards, which means that the ground state has changed 
from bonding to antibonding character. This ground 
state crossing is indicated by a green arrow in Fig.[TJa). 
The modulation of Ab a b with longitudinal magnetic fields 
is a consequence of the off-diagonal magnetic terms in 
Hi u tt ^ However, no such behavior is found in exper- 
iments, where Ab a b remains roughly constant with the 
field. This is shown by the symbols in Fig. [5J which rep- 
resent experimental data taken from Ref. [28|. Clearly, 
both Hmass (dashed-dotted line) and H K (solid line) suc- 
ceed in reproducing the approximately constant value of 

Abab- 



> 

s 




FIG. 3: Zeeman splitting as a function of the magnetic field. 
Dashed line: "Hiutt- Dashed-dotted line: Hmass- Solid line: 
H K . 

In order to discriminate between H mass and H K , in 
Figure El we compare the Zeeman splitting of the ground, 
A z , calculated with all three Hamiltonians and the ex- 
perimental values, represented by symbols. It can be seen 
that Hmass (dashed-dotted line) vastly overestimates the 
Zeeman splitting, while H K (solid line) offers the closest 
description. It is worth noting that the experimental val- 
ues of A z are even smaller than those predicted by H K . 
The inclusion of strain and piezoelectric effects may be 
relevant for a quantitatively improved description] 24 ' 42 

Last, we compare the envelope angular momentum ad- 
mixture obtained with the different Hamiltonians. In a 
four-band model, the hole states of cylindrical QDs are 
four-component spinors of the form: 



\F z ,n)= \M 

J z = -3/2,3/2 



(19) 



where n is the main quantum number and M = F z — J ' z 
is the envelope azimuthal angular momentum of a given 
component. We calculate the expectation value of the 
ground state envelope angular momentum at different 
fields and plot the results in Figure 0J At zero field 
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FIG. 4: Envelope angular momentum expectation value of 
the ground state as a function of the magnetic field. Dashed 
line: Hiutt- Dashed-dotted line: Hmass- Solid line: H K . 

(M) = 0.02, indicating that the ground state largest com- 
ponent is by far the HH (\J Z \ = 3/2) with M = 0, with 
a small admixture with finite M components. When the 
magnetic field is switched on, "Hiutt predicts a faster in- 
crease of (M) than rl maS s or rl K (the bump at B = 5 AT 
is due to the bonding-antibonding reversal) . Since the de- 
gree of envelope angular momentum admixture is critical 
in determining the effective g factor of holes and we 
have shown that only % K provides a consistent descrip- 
tion of the magnetic response, Fig. @] implies that the 
widely used Luttinger approximation is likely to overes- 
timate the g factor values in confined systems. 

IV. SUMMARY 

We have derived a multiband k-p Hamiltonian for va- 
lence holes confined in heterostructures subject to an ex- 



ternal magnetic field. The magnetic field has been imple- 
mented incorporating the minimal coupling p > p — qA 
prior to introduce the EFA. The inclusion of the remote 
bands has been considered through effective masses for 
the envelope function terms and effective g factors for 
the magnetic terms originating in the (unit cell) Bloch 
functions. For QDs, owing to the strong confinement, 
the latter have been replaced by the bare hole values, 
disregarding the influence of remote bands. 

The resulting Hamiltonian has been compared with 
the widely employed Luttinger approximation and our 
previous proposal. When tested against experimental 
data for InGaAs QDs under axial magnetic fields, the 
Hamiltonian presented in this work clearly outperforms 
the others. In particular, it succeeds in simultaneously 
describing the increase of the excitonic gap with the field, 
the constant splitting between bonding and antibonding 
states and the small Zeeman splitting observed in photo- 
luminescence experiments, with no fitting parameters. 

The Hamiltonian we have formulated is expected to 
improve current attempts to estimate the g factors of 
holes confined in QDs for spintronic, quantum informa- 
tion and optical applications. 
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